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Abstract 

\ We describe a method of constructing a sequence of phase coded waveforms with perfect autocorrelation 

in the presence of Doppler shift. The constituent waveforms are Golay complementary pairs which have perfect 
autocorrelation at zero Doppler but are sensitive to nonzero Doppler shifts. We extend this construction to multiple 
dimensions, in particular to radar polarimetry, where the two dimensions are realized by orthogonal polarizations. 
Here we determine a sequence of two-by-two Alamouti matrices where the entries involve Golay pairs and for 
which the sum of the matrix-valued ambiguity functions vanish at small Doppler shifts. The Prouhet-Thue-Morse 
sequence plays a key role in the construction of Doppler resilient sequences of Golay pairs. 

Index Terms 

Doppler resilient waveforms, Golay complementary sequences, perfect autocorrelation waveforms, Prouhet-Thue- 
Morse sequence, radar polarimetry. 

I. Introduction 

The value of perfect autocorrelation sequences in radar imaging is that their impulse-like autocorrelation function 
J> ■ can enable enhanced range resolution (e.g. see [l]-[6]). An important class of perfect autocorrelation sequences 
C*** . are complementary sequences introduced by Golay [6]. Golay complementary sequences have the property that the 
^ | sum of their autocorrelation functions vanishes at all (integer) delays other than zero. This means that the sum of 
£T) • the ambiguity functions (composite ambiguity function) of Golay complementary sequences is sidelobe free along 
O , the zero-Doppler axis, making them ideal for range imaging. 

The concept of complementary sequences was generalized to multiphase (or polyphase) sequences by Heimiller 
[7], Frank et al. [8]— [10], and Sivaswami [11], and to multiple complementary codes by Tseng and Liu [12]. Over 
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(j . the past five decades, the use of complementary and polyphase sequences (and related codes) have been widely 
£al | explored for radar imaging, e.g. see [1]— [18]. Recently, Deng [19] and Khan et al. [20] extended the use of polyphase 
sequences to orthogonal netted radar (a special case of MIMO radar), and Howard et al. [21] and Calderbank et 
al. [22] combined Golay complementary sequences with Alamouti signal processing to enable pulse compression 
C$ \ for multi-channel and fully polarimetric radar systems. Golay complementary sequences have also been advocated 
for the next generation guided radar (GUIDAR) systems [23]. 

Despite the attention they received from the radar engineering community, complementary and polyphase se- 
quences were somewhat ignored by communication engineers for many years, although their autocorrelation func- 
tions have as low sidelobes as the popular pseudo noise (PN) sequences. In fact, up until 1990, there were only 
a few articles on the use of complementary and polyphase sequences in communications, among which are the 
early work by Reed and Zetterberg [24] and the introduction of orthogonal complementary codes for synchronous 
spread spectrum multiuser communications by Suehiro and Hatori [25]. In 1990's, some researchers including 
Wilkinson and Jones [26], van Nee [27], and Ochiai and Imai [28] explored the use of Golay complementary 
sequences as codewords for OFDM, due to their small peak-to-mean envelope power ratio (PMEPR). However, 
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the major advances in this context are due to Davis and Jedwab [29] and Paterson [30], who derived tight bounds 
for the PMEPR of Golay complementary sequences and related codes from cosets of the generalized first-order 
Reed-Muller code. Construction of low PMEPR codes from cosets of the generalized first-order Reed-Muller code 
has also been considered by Schmidt [31] and Schmidt and Finger [32]. Complementary codes have also been 
employed as pilot signals for channel estimation in OFDM systems [33]. 

Recently, complementary and polyphase codes (in particular orthogonal complementary codes) have been advo- 
cated by Chen et al. [34], [35] and Tseng and Bell [36] for enabling interference-free (both multipath and multi- 
access) multicarrier CDMA. Other work in this context include the extension of complementary codes using the 
Zadoff-Chu sequence by Lu and Dubey [37] and cyclic shifted orthogonal complementary codes by Park and Jim 
[38]. In [39], orthogonal complementary codes have been used in the design of access-request packets for contention 
resolution in random-access wireless networks. 

Despite their many intriguing properties and recent theoretical advances, in practice a major barrier exists in 
adoption of complementary sequences for radar and communications; the perfect auto-correlation property of these 
sequences is extremely sensitive to Doppler shift. Although the shape of the composite ambiguity function of 
complementary sequences is ideal along the zero-Doppler axis, off the zero-Doppler axis it has large sidelobes 
in delay, which prevent unambiguous range imaging in radar or reliable detection in communications. Most 
generalizations of complementary sequences, including multiple complementary sequences and polyphase sequences 
suffer from the same problem to some degree. Examples of polyphase sequences that exhibit some tolerance to 
Doppler are Frank sequences [9], PI, P2, P3, and P4 sequences [18], PX sequences [40], and P(n, k) sequences 
[41], [42]. Sivaswami [43] has also proposed a class of near complementary codes, called subcomplementary codes, 
which exhibit some tolerance to Doppler shift. Subcomplementary codes consist of a set of N length- N sequences 
that are phase-modulated by a binary Hadamard matrix. The necessary and sufficient conditions for a set of phase- 
modulated sequences to be subcomplementary have been derived by Guey and Bell in [44]. The design of Doppler 
tolerant polyphase sequences has also been considered for MIMO radar. In [20], Khan et al. have used a harmonic 
phase structural constraint along with a numerical optimization method to design a set of polyphase sequences with 
resilience to Doppler shifts for orthogonal netted radar. 

In this paper, we present a novel and systematic way of designing a Doppler resilient sequence of Golay 
complementary waveforms for radar, for which the composite ambiguity function maintains ideal shape at small 
Doppler shifts. The idea is to determine a sequence of Golay pairs that annihilates the low-order terms of the Taylor 
expansion (around zero Doppler) of the composite ambiguity function. It turns out that the Prouhet-Thue-Morse 
sequence [45]-[48] plays a key role in determining the sequence of Doppler resilient Golay pairs. We then extend 
our analysis to the design of a Doppler resilient sequence of Alamouti waveform matrices of Golay pairs, for 
which the sum of the matrix-valued ambiguity functions vanishes at small Doppler shifts. Alamouti matrices of 
Golay waveforms have recently been shown [21], [22] to be useful for instantaneous radar polarimetry, which has 
the potential to significantly increase the performance of fully polarimetric radar systems, without increasing the 
receiver signal processing complexity beyond that of single channel matched filtering. Again, the Prouhet-Thue- 
Morse sequence plays a key role in determining the Doppler resilient sequence of Golay pairs. Finally, numerical 
examples are presented to demonstrate the perfect autocorrelation properties of Doppler resilient Golay pairs at 
small Doppler shifts. 

II. Golay Complementary Sequences 

Definition 1: Two length L unimodular sequences of complex numbers x[l] and y[l] are Golay complementary 
if the sum of their autocorrelation functions satisfies 



Then, x[l] and y[l] (or alternatively X(z) and Y(z)) are Golay complementary if X(z) and Y(z) satisfy 



coix k (x[l]) + corr fe (y[Z]) = 2L6 kfi , for k = —{L -!),■•■ , (L - 1), 



(1) 



where corr k {x[l\) is the autocorrelation of x[l] at lag k and 5 k fi is the Kronecker delta function. 
Let X(z) = Z{x[l]} and Y(z) = Z{y[l]} be the z-transforms of x[l] and y[l] so that 

X(z) = x[0] + x[l}z- 1 + . . . + x[L - l]*-^- 1 ) 
Y(z) = y[0] + y[l]z^ + ... + y[L- l]*"^" 1 ). 



(2) 



X(z)X(z) + Y(z)Y(z) = 2L 



(3) 
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or equivalently 

\\X(z)\\ 2 + \\Y{z)\\ 2 = 2L, (4) 

where X(z) = X*(l/z*) and Y(z) = Y*(l/z*) are the ^-transforms of x[l] = x*[-l] and y[l] = y*[-l], the time 
reversed complex conjugates of x[l] and y[l]. 

Henceforth we drop the discrete time index I from x[l] and y[l] and simply use x and y. We use the notation 
(a?, y) whenever x and y are Golay complementary and call (x, y) a Golay pair. From ([3]) it follows that if (x, y) 
is a Golay pair then (±x, ±y), (±x, ±y), and (±x, ±y) are also Golay pairs. 



A. Golay Pairs for Radar Detection 

Consider a single transmitter/single receiver radar system. Suppose Golay pairs (xq, xi), (x2, £3), ■ ■ ■ , (xn-2,xn-i) 
are transmitted over N pulse repetition intervals (PRIs) to interrogate a radar scene containing a stationary (relative 
to the transmitter and receiver) point target. Let R n (z) = Z{r n [k]} denote the z-transform of the radar return 
associated with the nth PRI. Then, the radar measurement equation can be written (in z-domain) as 

[R (z), Rn-^z)] = hz~ d ° [X (z), . . .,X N ^(z)} + [W (z), W N ^(z)} 

V V V (J ) 

r T (z) x T (z) w T (z) 

where the delay in z~ d ° depends on the target range vq and is given by do = [2ro/(cio)]> where to is the "chip" 
interval (time interval between two consecutive values in x[l] or y[l]), c is the speed of light, and [a] denotes the 
integer part of a. Without loss of generality, from hereon we assume that do = 0, centering the delay axis at the 
target location. The scalar h is the target scattering coefficient, which we assume to be proper complex normal 
with zero mean and variance 2af l , but fixed over the N RPIs. Elements of w T (z) are z-transforms of iid samples 
of proper complex white Gaussian noises with variance 2cr^. 

If we process the radar return vector r T (z) by a receiver vector of the form 

x(z) = [X (z),...,X JV -i(z)] T (6) 

then the receiver output will be 

U(z) = r T (z)x(z) = /ix T (z)x(z) + w T (z)x(z) 

(7) 

= NLh + w T (z)x(z), 
where the second equality follows by replacing x T ( z )x( 2 ) with 

x T (z)x(z) = (||Xo(z)|| 2 + WX^W 2 ) + ... + (\\X N - 2 (z)\\ 2 + \\X N ^(z)\\ 2 ) = NL. (8) 

The term *?{z)%{z) is the z-transform of the composite ambiguity function of Golay pairs {xq, x\), . . . , 
(xn-2,xn-i) along the zero-Doppler axis. We notice that x T (^)x(2;) is a constant, which means that the composite 
ambiguity function of (xo,xx), . . . , (xn-2, xn-i) vanishes at all (integer) delays along the zero-Doppler axis. 

Transforming {/J) back to the time domain, we have 

u[k] = Z- l {U{z)} = NLh5 kfi + n[k], (9) 

where n[k] is a proper complex white Gaussian noise with variance 2cr 2 = (A r L)2cr^. This shows that detecting a 
stationary point in range amounts to the following Gaussian hypothesis test 

f n ~ C7iV[0, 2a 2 n ] : H 

u[k] = { (10) 
[ NLh + n ~ CN[0, (2N 2 L 2 a 2 + 2a 2 n )} : Hi 

where CN[0, 2a 2 ] denotes the proper complex normal distribution with mean zero and variance 2a 2 . 
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Remark 1: In the above analysis, the radar return associated with each PRI is processed separately at the receiver, 
that is each radar return is correlated with its corresponding waveform and then all the correlator outputs are added 
together. Hence the receiver output (in time domain) is 

N-l 

u[k] = Yl xcorr fc (r n [£;'], x n [fc']) 

n=0 

N ~ l (11) 

= h Yl coiT k (x n [k']) +n[k] 

n=0 

= NLhS k ,o + n[k], 

where xcorr^ (r n [k'] , x n [k'\ ) is the cross-correlation between r n [k'] and x n [k'] at lag k. If we want to process all 
the PRIs together then we must correlate the augmented radar return r a [k], 

r a [k]=r [k]+r 1 [k-D] + ... + r N - 1 [k-(N-l)D] (12) 

with the augmented waveform x a [k], 

x a [k) = x [k] + Xl [k - D] + . . . + x N ^[k - {N - 1)D], (13) 

where D is the delay associated with a PRI. The receiver output in this case is 

u a [k] = xcorr fe (r a [k'\ , x a [k 1 ] ) 

N-l N-i N-l 

= h^2 corr k(x n [k']) + h ^ Y xc orr k (x nl [k' - n'D],x n »[k' — n"D}) + n a [k] 

n=0 n'=0n"=0 (\A\ 

n'^n" ' 

N-l N-l 

= NLh8 kfi + hY Yl xcorr fc (x n '[/c / - n'D],x n »[k' - n" D}) + n a [k], 

n'=0n"=0 
n'y^n" 

where n a [k] is a noise term. The cross terms xcorrj, (x n > [k' — n'D] , x n » [k' — n"D] ) result in range sidelobes whose 
peaks are offset by integer multiples of D from the origin k = 0. Thus, by processing each radar return separately 
as in (fTTl ) we can avoid range sidelobes caused by cross-correlations between different waveforms. However, the 
Doppler resolution will be limited by the time duration of a single waveform, whereas in the case where all the 
returns are processed together the Doppler resolution is enhanced due to having a longer transmit pulse. 



B. Effect of Doppler on Golay Pairs 

We now consider the case where the target moves at a constant speed, causing a Doppler shift of 9 [rad] between 
two consecutive PRIs. We assume that the radar PRI is short enough that during the N PRIs where the Golay pairs 
are transmitted the target range remains approximately the same. Then the composite radar measurement is given 

by 

r T {z, 6) = /ix T (z)D(0) + w T (z), (15) 
where D(0) is the following diagonal Doppler modulation matrix: 

D(0) = diag(l, e j9 , . . . , e j{N ~ l)e ). (16) 
If we now process the radar measurement vector r r (z, 9) using the receiver vector x(z) the receiver output will 

be 

U(z, 9) = r T (z, 9)x(z) = hG(z, 9) + w T (z)x(z), (17) 

where G(z, 9) is the z-transform of the composite ambiguity function of (xq, x±), . . . , (xn-2, <ejv-i), and is given 
by 

G(z, 9) = x T (z)D(0)x(z) = \\X (z)\\ 2 + e^WX^f + ... + e^ N -^ e \\X N ^(z)f . (18) 

We notice that off the zero-Doppler axis (9 ^ 0) the composite ambiguity function G(z, 9) is not sidelobe-free at 
integer delays. In fact, even small Doppler shifts can result in large sidelobes at integer delays. 
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One way to solve this problem is to use a bank of Doppler niters to estimate the unknown Doppler shift 9 and 
then compensate for the Doppler effect by post-multiplying (031 ) by T) H (9) (where H denotes Hermitian transpose) 
prior to applying x(z). However, since even a slight mismatch in Doppler can result in large sidelobes, we have to 
cover the possible Doppler range at a fine resolution, which requires the use of many Doppler filters. This motivates 
the question of whether it is possible to design Doppler resilient Golay pairs (xo,x\), . . . , (xn-2,xn~i) so that 
G(z,9) = J2n=o e " ,n6 'll^n(-2 ; )|| 2 ~ &z°, where a is constant, for a reasonable range of Doppler shifts 9. We are 
looking to construct the Golay pairs (xq, xi), . . . , {xn-2, xn-i) so that G(z, 6) (which is a two-sided polynomial 
of degree L — 1 in z^ 1 ) vanishes at every delay but zero. 



III. Doppler Resilient Golay Pairs 

In this section we consider the design of Doppler resilient sequences of Golay pairs. More precisely, we describe 
how to select Golay pairs (xo, xi), • • • , (xn-2, %n-i) so that in the Taylor expansion of G(z, 6) around 9 = the 
coefficients of all terms up to a certain order, say M, vanish at all nonzero delays. 

Consider the Taylor expansion of G(z, 9) around 9 = 0, i.e., 

oo 

G(z,9) = Y,C m (z)9 m , (19) 

where 

N-1 

C m {z) = Y,n m \\X n {z)\\ 2 , for m = 0,l,2,3,... (20) 

n=0 

In general, the coefficients C m {z), m = 1, 2, 3, . . . are two-sided polynomials in z^ 1 of the form 

L-l 

C m (z)= W** ™ = 1,2,3,... (21) 

l=-(L-l) 

For instance, the first coefficient C\ (z) is 

Ci(z) = 0||X (z)|| 2 + l||Xi(*)|| 2 + 2||X 2 (z)|| 2 + . .. + (N — l)\\X N ^{z)\\ 2 . (22) 

Noting that {xq,x\), . . . , (xn_2,xn-i) are Golay pairs we can simplify C\(z) as 

d(z) = N(N - 2)L/2 + WX^f + \\X 3 (z)\\ 2 + ... + \\X N ^(z)\\ 2 . (23) 

Each term of the form ||X 2 fc + i(z) || 2 = X 2 k+i(z)X2 k+1 (l/ z*) is a two-sided polynomial of degree L — 1 in the 
delay operator which can not be matched with any of the other terms, as we have already taken into account 
all the Golay pairs. Consequently, C\(z) is a two-sided polynomial in z~ x of the form C\{z) = Ylt=~(L-i) c l,l z ~ l - 
We wish to design the Golay pairs (xq, x\), . . . , (xn-2, xn-i) so that c\ \ vanish for all nonzero I. More generally, 
we wish to design (xo,x\), . . . , {xn-2, xn-\) so that in the Taylor expansion in ( fT9l the coefficients of all the 
terms up to a given order M vanish at all nonzero delays, i.e. c m ^ = 0, for all m (1 < m < M) and for all 
nonzero /. Although not necessary, we continue to carry the term O m ||Xo(z)|| 2 in writing C m (z) for reasons that 
will become clear. We note that there is no need to consider the zero-order term, as Cq(z) = NL. 



A. The Requirement that C\(z) Vanish at All Nonzero Delays 

To provide intuition, we first consider the case N = 2 2 = 4, where Golay pairs (xq,x\) and (x2, xs) are 
transmitted over four PRIs. Then, as the following calculation shows, C\(z) will vanish at all nonzero delays if 
the Golay pairs {xq,x\) and (x2,x 3 ) are selected such that [x\,x^) is also a Golay pair: 

C x {z) = 0||Xo(2)|| 2 + ||X 1 (z)|| 2 + 2||X 2 (z)|| 2 + 3||X 3 (z)|| 2 

lll^iWII 2 2x2L + l||X 3 (z)|| 2 (24) 

V v ' 

3 x 2L 
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The trick is to break 3 into 2 + 1, and then pair the extra ||X3(z)|| 2 with ||Xi(z)|| 2 . Note that it is easy to choose 
the pairs (xo,xi) and (x2,x 3 ) such that (x\,xs) is also a Golay pair. For example, let (x,y) be an arbitrary Golay 
pair, then (xo = x,xi = y), (x 2 = — y, X3 = x), and (xi = y,x 3 = x) are all Golay pairs. Other combinations of 
±x, ±x, ±y, and ±y are also possible. The calculation in (|24T > shows that it is possible to make C\(z) (M = 1) 
vanish at all nonzero delays with N = 2 1+1 Golay sequences xq, . . . , X3. 

B. The Requirement that C\{z) and (^2(2) Vanish at All Nonzero Delays 

It is easy to see that when ./V = 4 it is not possible to force C2(z) (M = 2) to zero at all nonzero delays. 
However, this is possible when N = 2 2+1 = 8. As the following calculations show, we can make both C\ (z) and 
62(2) vanish at all nonzero I if we select the Golay pairs (xo,xi), . . . , (xq,xj) such that (xi,xs), (X5, X7), and 
(x 3 ,x 7 ) are also Golay pairsQ 

Making C\{z) vanish: 

Cx(z) = 

0||X || 2 + lpi|| 2 + 2||X 2 || 2 + 3||X 3 || 2 + 4||X4|| 2 + 5||X 5 || 2 + 6||X 6 || 2 + 7||X 7 || 2 



2 x 2L + 4 x 2L + 6x 2L+ (25) 

7 



[(1-0) = 1]||X 1 || 2 [(3-2) = !] ||X 3 1| 2 [(5-4) = l]||X 5 || 2 [(7-6) = l]||^»-' 



3 x 2L 11 x 2L 

Making C2(z) vanish: 
C 2 (z) = 

2 ||X || 2 + l 2 ||Xi|| 2 + 2 2 ||X 2 || 2 + 3 2 ||X 3 || 2 + 4 2 ||X4|| 2 + 5 2 ||X 5 || 2 + 6 2 ||X 6 || 2 + 7 2 ||X 7 || 2 

4 x 2L + 16 x 2L + 36 x 2L+ 

[(1 2 -0 2 ) = 1]||X 1 || 2 [(3 2 - 2 2 ) = 5]||X 3 || 2 [(5 2 - 4 2 ) = 9]||X 5 || 2 [(7 2 - 6 2 ) = 13] ||X 7 || 2 (26) 

5 x 2L + 61 x 2L+ 

[(3 2 — 2 2 — l 2 + 2 ) =4]||X 3 || 2 [(7 2 - 6 2 - 5 2 + 4 2 ) = 4]||X 7 || 2 



[(0 2 + l 2 + 2 2 + . . . + 7 2 ) = 70] x 2L 

Note that it is easy to select the Golay pairs (xo,x{), . . . , (xq,x^) such that (x\,xs), (x5,x 7 ), and (x3,x 7 ) are 
also Golay pairs. For example, (xo = x, x\ = y), (x2 = —y,X3 = x), (X4 = — y, X5 = x), and (xe = x,xj = y), 
where (x, y) is an arbitrary Golay pair, satisfy all the extra Golay pair conditions. 

We notice that what allows us to make both C\(z) and C2(z) vanish at all nonzero I is the identity 

3"* _ 2 m - l m + m = 7 m - 6 m - 5 m + 4 m , for m = 1, 2 (27) 

or alternatively 

(0 m + 3 m + 5 m + 6 m )-(l m + 2 m + 4 m + 7 m ) = 0, for m = l,2, (28) 

where m = 1 and m = 2 correspond to the calculations for C\{z) and C^z), respectively. In other words, the 
reason C\{z) and C2(z) can be forced to zero at all nonzero delays is that the set S = {0,1,... ,7} can be 
partitioned into two disjoint subsets So = {0,3,5,6} and §1 = {1,2,4,7} whose elements satisfy d28l) . This is a 
special case of the Prouhet (or Prouhet-Tarry -Escort) problem [48], [49] which we will discuss in more detail later 
in this section. But for now we just note that So is the set of all numbers in S that correspond to the zeros in the 
length-8 Prouhet-Thue-Morse sequence (PTM) [45] -[48] 

(ajfc)Lo = 1 1 1 1, (29) 

and Si is the set of all numbers in S that correspond to the ones in (sfc)I=o- 

A key observation here is that the extra Golay pair conditions we had to introduce to make C\{z) and 62(2) 
vanish at all nonzero I are all associated with pairs of the form (x p ,x q ) where p and q are odd, and p G So and 
q € Si. This suggests a close connection between the Prouhet-Thue-Morse sequence and the way Golay sequences 
xq, x\, . . . , xjv_i must be paired. 



'in writing {25} and ( 126b we have dropped the argument 2 for simplicity. 
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C. The Requirement that C\{z) Through Cm{z) Vanish 

We now address the general problem of selecting the Golay pairs (xq, x\), . . . , (xn-2, sjv-i) to make C m (z), 
m = 1, 2, . . . , M vanish at all nonzero delays. We begin with some definitions and results related to the Prouhet- 
Thue-Morse sequence. 

Definition 2.[45]-[48] The Prouhet-Thue-Morse (PTM) sequence S = {s k )k>o over {0, 1} is defined by the 
following recursions: 

1) s = 

2) s 2k = Sk 

3) s 2 k+i = s k = l-s k 

for all k > 0, where s = 1 — s denotes the binary complement of s G {0, 1}. 
For example, the PTM sequence of length 32 is 

S = (s k )kL = 0110100110010110100101100110100 1. (30) 

Prouhet's prob!em.[48],[49] Let § = {0, 1, . . . ,N - 1} be the set of all integers between and N - 1. The 
Prouhet's problem (or Prouhet-Tarry-Escott problem) is the following. Given M, is it possible to partition S into 
two disjoint subsets So and Si such that ^ p m = q rn for all < m < M? Prouhet's proved that this is 

possible when N = 2 M+1 and that the partitions are identified by the PTM sequence. 
Theorem 1 (Prouhet).[48],[49] Let S = (s k )k>o be the PTM sequence. Define 

So ={p€S = {0,l,2,...,2 M+1 -l}| s p = 0} 

(31) 

Sx ={(?eS = {0,l,2,...,2 M+1 -l}| s q = l} 
Then, for any m with < m < M we have 

5> m =£<7"\ (32) 

Lemma 1. Let (xq,xi), . . . , (xn_2, x n-i)> N = 2 A/+1 be Golay pairs. Let X = {x p \p G So} and Xi = {x q \q G 
S q }. Then, neither Xo nor Xi contains any of the Golay pairs {xq : x\), . . . , {xn-2,xn-i)- 

Proof: The Golay pairs (xq, xi), . . . , (xn-2,xn-i) are of the form (x 2 k, ^2fc+i)> where k = 0, 1, . . . , N/2 — 1. 
From the definition of the PTM sequence we have s 2 k+i = Sk = s 2 k- Therefore, £ 2 fc an d a^fc+i cannot be in the 
same set.D 

Lemma 2. Assume that the Golay pairs (xo,x\), . . ., (xjy_ 2 , xjf-i), N = 2 AI+1 are such that all pairs of the 
form (x2k'+i € Xo,X2fe"+i € Xi), i.e., all pairs of the form (x 2 k'+i, ^2fc"+l) with 2A;' + 1 € So and 2k" + 1 G Si, 
are also Golay complementary. Then, 

\\X p {z)\\ 2 = \\X p ,(z)f and \\X q (z)\\ 2 = \\X q ,(z)\\ 2 (33) 

for all p,p' G So (i.e. for all x p , x p > G Xo) and for all q, q' G Si (i.e. for all x q , x q > G Xi), and all pairs of the form 
(x p G Xo,x g G Xi), i.e. all pairs (x p ,x q ) with p G So and q G Si, are Golay complementary. 

Proof: Assume p = 2k is even and p G So- Then q = 2k + 1 is odd and q G Si. We know that the pair 
(xp=2k £ Xo,Xq = 2fc+i G Xi) is Golay complementary, as all the original Golay pairs (xq,xi), . . . , (xn-2, %n-i) 
are of the form (x 2 k, x 2 k+i), hence 

\\X p (z)\\ 2 + \\X q (z)f = 2L. (34) 

Let p' G So and assume p' is odd. Then, since q = 2k + l G Si and all pairs of the form (x 2 k>+i G Xo, X2fc"+i G Xi) 
are Golay complementary (from our assumption), we have 

\\X p ,(z)\\ 2 + \\X q (z)\\ 2 = 2L. (35) 

Subtracting d33T ) from d34l ) gives 

\\X p (z)\\ 2 = \\X p ,(z)\\ 2 . (36) 

Since d36l ) is true for any even p G So and any odd p' G So it must be true for any p,p' G So, or equivalently any 
x p ,x p > G Xq. Similarly, we can prove that ||X g (z)|| 2 = ||X g '(z)|| 2 for all x q ,x q > G Xi. Since at least one element 



PEZESHKI et al: DOPPLER RESILIENT WAVEFORMS WITH PERFECT AUTOCORRELATION 



8 



from Xq forms a pair with one element in Xi (e.g. (xq,xi)) then all pairs of the form (x p £ Xo,x g £ Xi) must 
be Golay complementary. □ 

Remark 2: We note that to construct Golay pairs (xq,xi), . . . ,(xn-2,xn-i), N = 2 M+l that satisfy the 
conditions of Lemma 2 we can consider an arbitrary Golay pair (a;, y) and then arbitrarily choose x p £ Xo 
from the set {x, —x, x, —x} and x q £ Xi from the set {y, —y, y, —y}, for any p £ So an d any q £ Si. 

We now present the main result of this section by stating the following theorem. 

Theorem 2. The coefficients C±(z), . . . , Cm{z) in the Taylor expansion ( fl9l ) will vanish at all nonzero delays 
if the Golay pairs {xq,x\) : . . . , (xat_2, xn-\), N = 2 M+l are selected such that all pairs (x p ,x q ) where p and q 
are odd and p £ So an d q £ Si are also Golay complementary. 

Proof: From Lemma 2, we have ||X p (z)|| 2 = \\X p ,(z)\\ 2 for all p,p' £ S and ||X 9 (,z)|| 2 = \\X q ,(z)\\ 2 for all 
q, q' £ Si. Therefore, we can write C m {z) (1 < m < M) as 

C7 m (z) = £ n m ||X n (z)|| 2 = (J2 P m )\\X (z)\\ 2 + 9 m )ll^iWH 2 - (3V) 

n=o peSo geSi 

From the Prouhet theorem (Theorem 1), we have Yl P m = E ?™ = A where j3 is constant. Therefore, we have 

C m (z) = l3(\\Xo(z)f + ||Xi(z)|| 2 ) = 20L. (38) 



IV. Doppler Resilient Golay Pairs for Fully Polarimetric Radar Systems 

Fully polarimetric radar systems are capable of simultaneously transmitting and receiving on two orthogonal 
polarizations. The use of two orthogonal polarizations increases the degrees of freedom and can result in significant 
improvement in detection performance. Recently, Howard et al. [21] (also see [22]) proposed a novel approach to 
radar polarimetry that uses orthogonal polarization modes to provide essentially independent channels for viewing a 
target, and achieve diversity gain. Unlike conventional radar polarimetry, where polarized waveforms are transmitted 
sequentially and processed non-coherently, the approach in [21] allows for instantaneous radar polarimetry, where 
polarization modes are combined coherently on a pulse by pulse basis. Instantaneous radar polarimetry enables 
detection based on full polarimetric properties of the target and hence can provide better discrimination against 
clutter. When compared to a radar system with a singly-polarized transmitter and a singly-polarized receiver the 
instantaneous radar polarimetry can achieve the same detection performance (same false alarm and detection 
probabilities) with a substantially smaller transmit energy, or alternatively it can detect at substantially greater 
ranges for a given transmit energy [21]. 

A key ingredient of the approach in [21] is a unitary Alamouti matrix of Golay waveforms that has a perfect 
matrix-valued ambiguity function along the zero-Doppler axis. The unitary property of the waveform matrix allows 
for detection in range based on the full polarimetric properties of the target, without increasing the receiver signal 
processing complexity beyond that of single channel matched filtering. We show in this section that it is possible 
to design a sequence of Alamouti matrices of Golay waveforms, for which the sum of the matrix-valued ambiguity 
functions vanishes at all nonzero (integer) delays for small Doppler shifts. 

Figure \T\ shows the scattering model of the fully polarimetric radar system considered in [21], where hvH 
denotes the scattering coefficient into the vertical polarization channel from a horizontally polarized incident 
field. Howard et al. employ Alamouti signal processing [50] to coordinate the transmission of {N/2) Golay pairs 
(xo,xi), . . . , (xjv_2, xn-i) over vertical and horizontal polarizations during N PRIs. The waveform matrix is of 
the form 

{) [X^z) X (z) ... X 2k+1 (z) X 2k (z) ... X N ^(z) X N ^ 2 (z))' 

where different rows in X(z) correspond to vertical and horizontal polarizations, and different columns correspond 
to different time slots (PRIs). 

The radar measurement matrix R(z) for this transmission scheme is given by 



R(z) = HX(z)D(0) + WO), 



(40) 
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Fig. 1. Scattering model for a fully polarimetric radar system, with a dually-polarized transmit and a dually-polarized receive antenna. 



where H is the 2 by 2 target scattering matrix, with entries hyy, hyn, hnv> and hun, W(z) is a 2 by N noise 
matrix with entries that are iid proper complex normal with zero mean and variance 2<r^, and D(0) is the diagonal 
Doppler modulation matrix introduced in (fT6l ). 

If we process R(z) with a receiver matrix X(z) of the form 

'x Q {z) -X x {z) ... X 2k (z) -X 2k+l (z) ... X N - 2 {z) -X N ^(z)\ T 
Xx{z) X (z) ... X 2k+1 (z) X 2k (z) ... X N ^(z) X N ^ 2 (z) J 



X(z) 

then the receiver output will be 



(41) 



U(z, 9) = K(z)X(z) = HG(z, 0) + W(*)X(z) 



(42) 



where the term G(z, 0) = X(z)D(0)X(z) can be viewed as the z-transform of a matrix-valued ambiguity function 
for X(z). Along the zero-Doppler axis, where D(0 = 0) = I, due to the interplay between Alamouti signal 
processing and the Golay property, the term G(z, 9) reduces to 



G(z,0) = X{z)X{z) 



/N-l 

' £ \\X n (z)\\ 2 = NL 

n=0 



N/2-1 _ \ 

Z (l-l)X 2k (z)X 2k+1 (z)=0 X 

k=0 



N/2-1 

£ (l-l)X 2k (z)X 2k+1 (z)=0 

k=0 



N-l 

£ \\Xn(z) 
n=0 



NL 



(43) 



This shows that X(z) has a perfect matrix-valued ambiguity function along the zero-Doppler axis; that is along the 
zero-Doppler axis G(z, 9) vanishes at all nonzero (integer) delays, and is unitary at zero-delay. A consequence of 
d43l ) is that detecting a point target in range reduces to a simple Gaussian hypothesis test, for which the likelihood 
ratio detector is the same as an energy detector. However, off the zero-Doppler axis the property in (l43l) no longer 
holds, and the elements of the matrix-valued ambiguity function G(z, 9) can have large sidelobes, even at small 
Doppler shifts. 

We consider how the Golay pairs {xq,x\), . . . , (xn-2, %n-i) must be selected so that for small Doppler shifts 
we have 

'(?i(z,0) G 2 {z,9)\ _ (NL 
y G 2 (z,9) G 1 (z,9)J ~{0 NL, 



G(z,0) = X(z)D(0)X(z) 



(44) 



where 



N-l 



G x {z,9) = £ e^ e \\X n {z) 

n=0 



(45) 



\X (z)\\ 2 + e^WX^f + ... + e^ N -^\\X N ^(z)\\ 2 



and 



N/2-1 



G 2 (z,9)= £ (e^ ke - e^ 2k +^)X 2k (z)X 2k+1 (z) 



k=0 



(46) 



(1 - e> e )X (z)X 1 (z) + ... + (e^-W - e^ N -^ )X N - 2 (z)X N ^(z). 
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The diagonal term of G(z, 9), i.e., G\(z, 9), is equal to the single channel composite ambiguity function G(z, 9) 
in (fT8l) . Therefore, we can use Theorem 2 to design the Golay pairs {xq,x\), . . . , (xat_2, xn-i), N = 2 M+1 such 
that in the Taylor expansion ( f!9l ) the coefficients C m (z), m = 1, 2, . . . , M vanish at all nonzero delays. Thus, from 
now on we only discuss how the off-diagonal term G 2 (z, 9) can be forced to zero for small Doppler shifts. 

Consider the Taylor expansion of G 2 (z, 9) around 9 = 0, i.e., 

00 

G 2 (z,9) = B m (z)9 m , (47) 

m=0 

where 

N/2-1 

B m {z) = m) m -(2k + l) m )X 2k (z)X 2k+1 (z), for m = 0,1,2,.... (48) 

k=0 

In general, the coefficients B m {z), m = 1, 2, 3, . . . are two-sided polynomials in z" 1 of the form 

L-l 

B m {z)= Yl h m^~\ rn= 1,2,3,... (49) 

1=-{L-X) 

For instance, the first coefficient B\ (z) is 

Bx{z) = (0 - l)X (*)*iCz) + . . . + {{N — 2) — [N — \))X N - 2 {z)X N .. 1 (z). (50) 

Each term of the form X 2 k{z)X 2 j t+ i(z) in (|50l ) is a two-sided polynomial of degree L — 1 in and since in 
general the terms X 2 k(z)X 2 k + i{z) for different values of k do not cancel each other, B\(z) is also a two-sided 
polynomial of degree L — 1 in z^ 1 . 

Suppose that the Golay pairs (xq,xi), . . . ,(xn- 2 ,xn-i), N = 2 M+l satisfy the conditions of Theorem 2 
so that Ci(z),... ,Cm(z) vanish at all nonzero delays. We wish to determine the extra conditions required for 
(xq, xi), . . . , (xn-2, %n-i) to force B\{z), . . . , Bm{z) to zero at all delays. As we show, again the PTM sequence 
is the key to finding the zero-forcing conditions. The zero-order term Bq(z) is always zero and hence we do not 
consider it in our discussion. 



A. The Requirement that B\{z) Vanish 

Again, to gain intuition, we first consider the case N = 2 2 = 4. Then, as the following calculation shows, B±(z) 
will vanish if the Golay pairs (xq,xi) and (^2,^3) are selected so that Xq(z)X\(z) = —X 2 {z)X^{z): 

B x {z) = (0 - l)X (z)Xx(z) + (2 - 3)X 2 (z)X 3 (z) 

^ v ' 

-(2-3)Xo(z)X 1 (z) (51) 





In summary, to make C\{z) vanish at all nonzero delays and to force B\{z) to zero at the same time, the Golay 
pairs (xq,x\) and (x 2 , £3) must be selected such that (x\,Xz) is also a Golay pair and Xq(z)Xi(z) = —X 2 (z)X^{z). 
If we let (x,y) be an arbitrary Golay pair then it is easy to see that (xq = x,x\ = y), (x 2 = —y,x 3 = x) satisfy 
these conditions (other choices are also possible). The Alamouti waveform matrix X(z) for this choice of Golay 
pairs is given by 

Xh _ r . ----- -X 1 (z) = -Y(z) X 2 (z) = -Y(z) -X 3 (z)= 

X 1 " ' X (z) = X(z) X 3 (z) = X(z) X 2 (z)= r 
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B. The Requirement that B\{z) and B 2 (z) Vanish: 

Let us now consider the case N = 2 3 = 8. Then, as the following calculations show, both B\(z) and B2{z) will 
vanish if we select (xq,x\), . . . , (xq,x 7 ) such that Xq(z)X\(z) = — X 2 (z)X 3 (z) = — X±{z)X*,(z) = Xq(z)X 7 (z). 
Making B\(z) vanish: 

B ^ = 

(0-l)X X 1 + (2-3)X 2 X 3 + (4-5)X 4 X 5 + (6 - 7)X 6 X 7 



[(0 - 1) = -1]X X 1 [-{2-3) = l]X X 1 [-(4-5) = 1}X X 1 [(6 - 7) = -l]X X 1 

N v 



Making B2(z) vanish: 

{0 2 -l 2 )X X 1 + (2 2 - 3 2 )X 2 X 3 + (4 2 - 5 2 )X 4 X 5 + (6 2 - 7 2 )X 6 X 7 



(53) 



[(0 2 — l 2 ) = —1]XqXi [-(2 2 - 3 2 ) = 5}X X 1 [-(4 2 - 5 2 ) = 9}X X 1 [(6 2 - 7 2 ) = -1S]X Q X 1 (54 ) 



[0 2 - l 2 - 2 2 + 3 2 = 4]X Xi [-4 2 + 5 2 + 6 2 - 7 2 = -4]X Xi 





In writing d53l ) and (|54l ) we have dropped the argument z for simplicity. 

In summary, to make C\{z) and C%(z) vanish at all nonzero delays and to force B\{z) and B 2 {z) to zero at 
the same time, the Golay pairs [xq, x\), . . . , (xq,x 7 ) must satisfy the conditions of Theorem 2, and the within-pair 
cross-spectral densities must satisfy 

X (z)Xi(z) = -X 2 (z)X 3 {z) = -X 4 (z)X 5 (z) = X 6 (z)X 7 (z). (55) 

It is easy to see that the Golay pairs in the following waveform matrix X(z) satisfy all these conditions: 

X(z) = (Xo(z) Xi(z) X x (z) X (z)) , (56) 

where Xo(z) and Xi(z) are given by 

Xo(2) = (m -?i>i\ md Xl(z) = ( -jw -^A _ (57) 

X(z) ) \X(z) -Y{z))' 

and (x, y) is an arbitrary Golay pair. 

The trick in forcing B\ (z) and B 2 (z) to zero is to cleverly select the signs of the cross-correlation functions 
(cross-spectral densities) between the two sequences in every Golay pair relative to the cross-correlation function 
(cross-spectral density) for xq and x\. If we let and 1 correspond to the positive and negative signs respectively, 
we observe that the sequence of signs in (|55l ) corresponds to the length-4 PTM sequence. In the next section, we 
show that the PTM sequence is in fact the right sequence for specifying the relative signs of the cross-correlation 
functions between the Golay sequences in each Golay pair. 

Remark 3: Representing Xq(z) and Xi(z) by and 1 respectively, we notice that the placements of Xo(-z) and 
Xi(z) in X(z) are also determined by the length-4 PTM sequence. 



C. The Requirement that B\(z) Through Bm(z) Vanish 

We now consider the general case N = 2 M+1 where Golay pairs (xq, x\), . . ., (xn-2, xn-i) are used to construct 
a Doppler resilient waveform matrix X(z). We have the following theorem. 

Theorem 3: Let N = 2 AI+1 and let (xq,xi), . . ., (xn-2, %n-i) be Golay pairs. Then, for any m between 1 
and M, B m (z) will vanish at all delays if 

X 2k {z)X 2k+1 {z) = (-IY*Xq{z)X x {z), for all k, < k < N/2 - 1, (58) 

where is the kth element in the PTM sequence. 
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Proof: For any m (1 < m < M), B m (z) may be written as 



N/2-1 

B m (z) = E ((2*0 m " (2* + l) m )X 2fe (z)X 2fe+1 (z) 

fc=0 



N/2-1 

£ (-l) s *((2£;) m - (2A; + l) m ) 

fc=0 



(59) 



X (z)*i(z), 



where the second equality in (|59l follows by replacing X2fc(z)X2fc+i(z) with (— l) Sk Xq{z)X\{z). Since in the 
PTM sequence Sk = S2k = S2fc+i« we can rewrite d59l as 



B m (z) 



N/2-1 

E (-l) S2fc (2fc) m - (-l) S2fc + 1 (2fc + l) r 

fc=0 
jV/2-l 

E (-l) S2fc (2A;) m + (-l) s ^+i(2fc + l) f 

fc=0 



XnfzlXi (z) 



X {z)X 1 {z) 



(60) 



2V-1 

E (-1)'**;" 



X (z)Xi(z). 



jV-1 



However, from the Prouhet theorem (Theorem 1), it is easy to see that ^ (— l) Sk k m = 0, and therefore B m (z) = 

on 

Finally, we note that it is always possible to find Golay pairs (xq, x\), . . . , (xn-2, xn-i) that satisfy the conditions 
of both Theorem 2 and Theorem 3. Suppose (xo, xi), . . . , (xat_2, £v-i) are built from an arbitrary Golay pair (x, y) 
(as explained in Section ITITb to satisfy the conditions of Theorem 2. Then, we can apply the time reverse operator 
and change the sign of the elements within the pairs to satisfy the conditions of Theorem 3, as the Golay property 
is invariant to time reversal and changes in the signs of the Golay sequences within a pair. 



V. Numerical Examples 

In this section, we present numerical examples to verify the results of Sections HIT] and [IV] and compare our 
Doppler resilient design to a conventional scheme, where the same Golay pair is repeated. 



A. Single Channel Radar System 

We first consider the case of a single channel radar system. In this case, the composite ambiguity function 
G(z, 6) is given by (|T8T > and has a Taylor expansion of the form ( fT9l . Following Theorem 2, we coordinate 
the transmission of eight Golay pairs (xq, x\), . . . , (sy, X15) over iV = 16 PRIs to make the Taylor expansion 
coefficients C\(z), . . . , C%(z) (M = 3) vanish at all nonzero delays. Starting from a Golay pair (x, y), it is easy to 
verify that the eight Golay pairs in the following waveform vector x T (z) satisfy the conditions of Theorem 2: 

x T (z) = (xJ(z) xf(z) xf(z) xj(z) xj(z) xj(z) *g(z) (61) 

where xj(z) = [X(z) Y(z)] and xf (z) = [-Y(z) X(z)}. 

Remark 4: Representing xj(z) and xf(z) by and 1 respectively, we notice that the placements of Xo (z) and 
xf (z) in x T (z) are determined by the length-8 PTM sequence. 

We compare the Doppler resilient transmission scheme in doTT) with a conventional transmission scheme, where 
the same Golay pair (xq = x,x± = y) is transmitted during all PRIs, resulting in a waveform vector xj(z) of the 
form 

xj(z) = (xj(z) xj(z) xj(z) xj(z) xj(z) xj(z) xj(z) xj(z)), (62) 
with the composite ambiguity function 

G c (z,9) = xJ(z)D(0)x c (z) 

(63) 

= \\X (z)\\ 2 + e^HX^z)!! 2 + . . . + e^- 2 ) e ||X (z)|| 2 + e^'^X^f . 
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The pair (x,y) used in constructing x T \z) and xj(z) can be any Golay pair. Here, we choose (x,y) to be the 



following length-8 (L = 8) Golay pair: 

x[l] ={1,1,-1,1,1,1,1,-1} 
y[l] ={-1,-1,1,-1,1,1,1,-1} 



X(z) 
Y(z) 



l + z~ l 
-1 - z~ 



z~ 2 + z~ 3 + z~ A 



z~ 3 + z~ 4 + z~ 5 + z~ 6 



y-7 



(64) 



Referring to the Taylor expansion of G(z, 9) in ( fl9l ), the coefficients C\{z), 6*2(2;), and Cs(z) are each two-sided 
polynomials of degree L — 1 = 7 in z^ 1 of the form (f2TT ). Figures I3a)-(c) show the plots of the magnitudes of 
the coefficients c mj i, m = 1, 2, 3 of these polynomials versus delay index I. The plots show that C\(z), C2(z), and 
Cs(z) indeed vanish at all nonzero delays. 

Figures [Ha),(b) show the plots of the composite ambiguity functions G(z,9) and G c (z,9) versus delay index 
I and Doppler shift 8. Comparison of G(z,9) and G c {z,9) at Doppler shifts 9 = 0.025 rad, 8 = 0.05 rad, and 
9 = 0.075 rad is provided in Figs. IUa)-(c), where the solid lines correspond to G(z, 9) (Doppler resilient scheme) 
and the dashed lines correspond to G c (z, 9) (conventional scheme). We notice that the peaks of the range sidelobes 
of G{z, 9) are at least 24 dB (for 9 = 0.025 rad), 28 dB (for 9 = 0.05 rad), and 29 dB (for 9 = 0.075 rad) smaller 
than those of G c (z,9). These plots clearly show the Doppler resilience of the waveform vector in (|6Tb . 

Remark 5: For a radar with carrier frequency /o = 2.5 GHz and PRI= 100 fisec, the Doppler shift range of to 
0.05 rad (0.075 rad) corresponds to a maximum target speed of V ~ 35 kmph (50 kmph). To cover a larger speed 
range we can use our design with a bank of Doppler filters to provide Doppler resilience within an interval around 
the Doppler frequency associated with each filter. 
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(a) h;| 



(b) |c 2 ,| 



(c) |c 3 ;| 



Fig. 2. The plots of the magnitudes of the coefficients c m> ! of two-sided polynomials C m (z) = c mjz l , m = 1,2,3 versus 

l=-{L-\) 

delay index /: (a) jci^l, (b) \c 2 ,i\, and (c) |c 3 j|. 




(a)G(z,9) (b)G c (z,9) 

Fig. 3. (a) The plot of the composite ambiguity function G(z, 9) (corresponding to the Doppler resilient transmission scheme) versus 
delay index I and Doppler shift 6, (b) the plot of the composite ambiguity function G c (z,0) (corresponding to the conventional transmission 
scheme) versus delay index I and Doppler shift 9. 
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Fig. 4. Comparison of the composite ambiguity functions G[z, 9) and G c (z, 9) at Doppler shifts (a) 9 — 0.025 rad, (b) 9 — 0.05 rad, and 
(c) 9 = 0.075 rad. 



B. Fully Polarimetric Radar System 

We now consider the matrix-valued composite ambiguity function G(z,9) in (l44l ). corresponding to the fully 
polarimetric radar system described in Section [IVl Following Theorems 2 and 3, we coordinate the transmission 
of eight Golay pairs (xq, Xi), . . . , (xu, X15) across vertical and horizontal polarizations and over N = 16 PRIs, so 
that in the Taylor expansions of G\{z, 9) (the diagonal element of G(z, 8)) and G 2 (z, 9) (the off-diagonal element 
of G(z,9)) the coefficients C\(z), C 2 (z), and C%(z) vanish at all nonzero delays and B\(z), B 2 {z), and B%{z) 
vanish at all delays. Letting ~Kq(z) and Xi(z) be the Alamouti matrices in (1571 ), then it is easy to check that the 
Golay pairs in the following waveform matrix X(z) satisfy all the conditions of Theorems 2 and 3: 

X(*) = (Xo(z) Xi(z) Xi(z) Xo(z) Xi(z) X (z) X (z) Xi(z)) . (65) 

Remark 6: Representing Xo(^) and Xi(z) by and 1 respectively, we notice that the placements of Xo(^) and 
Xi(z) in X(z) are determined by the length-8 PTM sequence. 

We compare the Doppler resilient transmission scheme in (|65T ) with a conventional transmission scheme, where 
the Alamouti waveform matrix built from a single Golay pair (xq = x, x\ = y) is repeated and the waveform 
matrix is of the form 

X c (z) = (Xq(z) Xo(z) Xo(z) X (2) Xq(z) Xq(z) X (z) X (z)) . 

The matrix- valued composite ambiguity function of X c (z) is given by 

G c (M)=XeWD(fl)X c ( Z ) = {c c2 {z,0) G cl (z,9))> 

where 

G cl (z, 6) = \\X {z)f + ^IIXjC^H 2 + . . . + e^ N - 2 ) d \\X {zW + e^"^^^!! 2 

and 

G c2 (z,9) = {l-ei e + ... + S N ~ 2 ) e - e^-^X^X^z). (69) 

The Golay pair (x,y) used in building both X(z) and X c (z) is the length-8 Golay pair in d64l ). 

We notice that the diagonal elements of G(z,9) and G c (z, 9), i.e., G\(z, 9) and G c \(z,9), are equal to the 
single channel composite ambiguity functions G(z,9) and G c (z,9), respectively. Therefore, the plots in Fig. |2] 
through Fig. [4] also apply for comparing G\{z,9) and G&(z, 9). Thus in this example, we only need to consider 
the off-diagonal elements of G(z,9) and G c (z,9), i.e., G 2 {z,9) and G c2 (z,9). 

Referring to the Taylor expansion of G 2 (z,9) in ( |47T ), the coefficients B\(z), B 2 {z), and B${z) are each two- 
sided polynomials of degree L — 1 = 7 in z^ 1 of the form (|49l ). Figures I3 a )-( c ) show the plots of the magnitudes 
of the coefficients b m i, m = 1,2,3 of these polynomials versus delay index I. We notice that B\{z), B 2 (z), and 
B3 (z) indeed vanish at all delays. 

Figures |Ma),(b) show the plots of the off-diagonal elements G 2 (z,9) and G c2 (z,9) versus delay index I and 
Doppler shift 9. Comparison of G 2 (z, 9) and G c2 (z, 9) at Doppler shifts 9 = 0.025 rad, 9 = 0.05 rad, and 9 = 0.075 



(66) 

(67) 
(68) 
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rad is provided in Figs. |7{a)-(c), where the solid lines correspond to G 2 (z,9) (Doppler resilient scheme) and the 
dashed lines correspond to G c2 (z,9) (conventional scheme). We notice that the peaks of the range sidelobes of 
G 2 {z, 9) are at least 24 dB (for 9 = 0.025 rad), 12 dB (for 9 = 0.05 rad), and 5 dB (for 9 = 0.075 rad) smaller than 
those of G c2 (z, 9). These plots together with the plots in Figs. Ua)-(c) (corresponding to the diagonal elements of 
G(z, 9) and G c (z,9)) show the Doppler resilience of the waveform matrix in 



: o ; a 



.5 o i S e 



^ 5 3 



(a) |6 



(b) 1 6 2 ,i I 



(c) |63,I I 



Fig. 5. The plots of the magnitudes of the coefficients b m .i of two-sided polynomials B m (z) = ^2 b m .iz , m = 1,2, 3 versus 

!=-(£-!) 

delay index I: (a) |6i,j|, (b) |& 2 ,i|, and (c) |6 3ii |. 




(a)G 2 (z,0) (b)G c2 (z,9) 

Fig. 6. (a) The plot of the off-diagonal element G2 (.£,#) (corresponding to the Doppler resilient transmission scheme) versus delay index 
I and Doppler shift 9, (b) the plot of the off-diagonal element G C 2(z,8) (corresponding to the conventional transmission scheme) versus 
delay index I and Doppler shift 6. 



VI. Conclusions 

We have constructed a Doppler resilient sequence of Golay complementary waveforms with perfect autocorrelation 
at small Doppler shifts, and extended our results to the design of Doppler resilient Alamouti waveform matrices of 
Golay pairs for instantaneous radar polarimetry. The main contribution is a method for selecting Golay complemen- 
tary sequences to force the low-order terms of the Taylor expansion of a composite ambiguity function (or Doppler 
modulated autocorrelation sum) to zero. The Prouhet-Thue-Morse sequence was found to be the key to selecting the 
Doppler resilient Golay pairs. Numerical examples were presented, demonstrating the perfect correlation properties 
of Doppler resilient Golay pairs at small Doppler shifts. 
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Fig. 7. Comparison of the off-diagonal elements G2(z,8) and G C 2(z,8) of the matrix- valued ambiguity functions for the conventional and 
Doppler resilient schemes at Doppler shifts (a) 6 — 0.025 rad, (b) 8 = 0.05 rad, and (c) 6 = 0.075 rad. 
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